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A lower bound is given for the quantity Xi/X„, and an upper bound for the quantity Xi — 
X„, where Xi and Xn are respectively the greatest and least characteristic roots of a matrix 
with positive roots. The bounds involve the first and second coefficients of the character- 
istic equation of the matrix. 



Suppose A= (aij) is a uonsiiio-ular nXn niatrix with 
characteristic roots X], Xo, . . ., ^n, so ordered that 
|Xil>|X2|> . . . >|X,|. The quantity |Xj|/|X,J pro- 
vides a rough measure of the probable error in the 
computation of the inverse of A; it has been called by 
J. Todd [1, 2, 3] ^ the P-condition nund)er of A and 
m^ay be denoted P{A). Von Neumann and (joldstine 
[4] have shown that if A is S3'mmetric and positive 
definite (in which case the X/ are all positive), then the 
error in the inverse of A computed by a certain elimi- 
nation method, can be bounded by a quantity pro- 
portional to P(A); if A is not symmetric positive 
definite, the error can be bounded by a quantity pro- 
portional to P(AA'). 

We shall restrict our consideration to matrices 
whose roots are all positive. For these, P. J. Davis, 
E. V. Haynsworth, and M. Marcus [5] obtained 
bounds on P involving det A and one other sym- 
metric function of the roots of A. If the character- 
istic polvnomial of A is p(x)=x''—C\x''~'^-{-C2X"'~^ + 
. . . +{-lYCr, and we set Di=(n''Cn)IC'i, they 
showed that 



A^ 



1 -^-i-vr=:F/ 



(1) 



and they found similar inequalities involving Cn 
(=det A) and any other one of the Ci. 

In many cases, however, Cn is not known and it is 
in general difficult to calculate. C](== trace A) is 
easy to calculate, and C2 can in general be calculated 
more easily than can Cn since it is the sum of 
[n(7i— 1)]/2 determinants of order two. In this paper 
we present (in Theorems 1 and 1') a lower bound for 
P in terms of Ci and C2; an attempt to obtain a cor- 
responding upper boimd fails, but leads to an inequal- 
ity (Theorem 2) on Xi— X„, the ^^spread'' of the roots 
of A. Finally we apply the method of proof of The- 
orem 1 to obtain an improvement of the low^er bound 
in (1). 



' Figures in brackets indicate the literature references at the end of this paper. 



If we set Ck={k)Sk(\u X., . • ., X,0, then SU"" is 
the K'ih synnnetric mean of Xi, . . ., X„. The Sk 
satisfy [6] the incnpialities 



S,>SV'>SV'> . . . >S'J\ 



(2) 



Setting tJLK=v^^ we have i^=Mi^M2^M3^ • • • 

A// 



and 



Let 



Soi^^l, M2, . . ., Mn) ^^'.(Xi, . . ., X„) 



Il{X2, X^, . . ., Xn-\) 



S2{P, X2,Xs, . . ,,Xn-l,l) 



and 



/i(P)=: max n(x2, ^3, . . . , ^..-1) 

1<2-,<P 



J2{P)= inin R{X2, . . . , ^„_i). 

1<3:,<P 



Then/i and/2 can be seen to be increasing in P, and 

(3) 



/2(P)<f</,(P). 



Thus the right half of (3) should provide a lower 
bound for P, while the left half should provide an 
upper bound. We first calculate y2(P) '- 

By direct calculation we can show that (d^P)/da:f 
is iionnegative at all points for each i. Therefore R 
attains its maximum at a point where each Xi is 
either 1 or P. Letting Rk{P) denote the value of R 
when K—l of X2^ x^, . . ., Xn-\ are equal to P and the 
remainder are equal to one, we find that Rk is equal 
to: 

n-l (KP+n-Ky . 

n (KP+n-Ky-iKP' + n-K)' ^ ^ 
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This rational function of K attains its maximum at 
n 



K- 



"P+l 



and so we obtain 



S\ n — 
S2~ n 



which is equivalent to: 
Theorem 1: 



1 



(p+iy 



(5) 



4Pn 



l+^Jl-A 



<P, where A= 



^-c^ (^-1) 



SI 



The upper bound (5) can be sharpened, since in 
fact we need only consider integer values of 
Kin (4). It can be shown that i?i=niax Hk if P^ 

>[(n-l)(n-2)]/2. Thus if P'>[{n-l){n-2)]/2, 
SlfS2<Bi(P)] and so, setting p=Si/S2, we have: 



[(p-l)+Vp(p-l)]n + l<P. 



(6) 



Now if P2<[(/i-l)(7i-2)]/2, then p<niax Rk{P) 

K 

=Bk^{P), say; and since each Rk(P) is a strictly 
increasing function of P, p</?^*([(7i— l)(n— 2)]/2) 



</?i([(7i-l)(7i-2)]/2)andso 



V 



{n—l){n — 2) 



>Ri-Kp)^ 



However it can easily be seen that if p>(V2+|)/ 
(,/^+l), then 7?.-(p)>-^EIMEI); 



and so we 



// p>(V2+f)/(V2+l), then [(p-1) 



conclude that: 

Theorem V: 
+Vp(p-l)]n+l<P. 

This lower bound is better than the previous one, 
whenever it applies. 

The attempt to derive an upper bound for P 
from the left half of (3) fails, because/2 (P) approaches 
1 uniformly in P as n increases, and so the inequality 
f2iP) ^81/82 will in most cases hold for all values of 
P. We can, however, by considering the minimum 

of SK\ y2, Vs, . . ., Vn^U K)—S2{K y2, . . ., Vn-h 

\n) subject to the condition Xi>?/2>2/3^. • -^Vn-i 
>\ny obtain an upper bound for the spread of the 
roots of A. Calling the above function of 7/2, . . ., 
yn-i D, we evaluate its minimum directly by observ- 
ing that c)D/dyi=-2/[n(n-l)]{S,-yi) and d'D/ 
^yi^y3=^/n^ if i=j and —2/[n\n—l)] if i^^j; 
hj=2, 3, . . ., n—l. The (n—2)X(n-2) matrix 
(dij), with dij=d^D/c>y,dyj is symmetric, and by 
Gerschgorin's theorem each of its eigenvalues lies in 



the circle \Z—2/n^\<2/n' 



(m 



and so is positive. 



Thus (dij) is positive definite, and, setting each 
c>D/dyi equal to zero, we find that D is at a minimum 
when each ?/, is equal to {\i-]-\n)/2, and that the 
minimal value is [(Xi— Xw)^]/[27i(7i— 1)]. We may 

then conclude: 

Theorem 2: Xi-X, <V2/z^(7i-l)(S?— &). 



By the method of Theorem 1, it is possible to 
sharpen the left side of inequality (1). Following 
the notation of [5] we write 



Ul{Xi,X2, . . ., ^;^) = 






/ X,+ . . . +Xn \ ^ 



and seek an upper bound for Di subject to the 
condition P=Xi>X2> . . . >:Cw=l. As in the proof 
of Theorem 1, we show that d^M/5x?>0 for i==2, 3, 
. . ., n—l and finally obtain the relation 






(7) 



which leads to the inequality: 
Theorem 3: 



< 



P-1 



logP 



i logP * 

This inequality yields a lower bound for P which 
is always higher than that given by (1) (as may be 
seen by comparing the proof of Theorem 3 with the 
proof of (1) in [5]). However it is cumbersome. 
It can be simplified (and somewhat weakened) as 
follows: Since xe^^^=a, a>6, implies that x>a—e-{-l, 

j_ 
and e/Di''>e by (2), we may conclude that (P— 1)/ 

]_ 
logP>6/Z)i^ — 6 + 1. Denoting this last quantity by 
A, we have (P— l)/log P>A, which has as an im- 
mediate consequence: 

j_ 

Theorem 3': P>Alog A-l] A=e/D,'^-e+l, 
This lower bound is most often, though not always, 
better than that given in (1). 
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